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ABSTRACT 

'n— ' . In the present work, we investigate the tension in the recent Type la supernovae (SNIa) datasets 

Constitution and Union. We show that they are in tension not only with the observations of the 
' cosmic microwave background (CMB) anisotropy and the baryon acoustic oscillations (BAO), but 

also with other SNIa datasets such as Davis and SNLS. Then, we find the main sources responsible 
for the tension. Further, we make this more robust by employing the method of random truncation. 
Based on the results of this work, we suggest two truncated versions of the Union and Constitution 
datasets, namely the UnionT and ConstitutionT SNIa samples, whose behaviors are more regular. 
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I. INTRODUCTION 

The current accelerated expansion of our universe (see e.g. [l|) has been one of the most active fields 
in modern cosmology since its discovery from the observations of Type la supernovae (SNIa) [2] . Later 
on, the observations of the cosmic microwave background (CMB) anisotropy 13] and the large-scale 
structure (LSS) Q confirmed this discovery. Although there are already many observational methods to 
date, SNIa have been proven to be one of the most powerful tools to probe this mysterious phenomenon. 
In the passed decade, many SNIa datasets have been released, while the number and quality of SNIa 
have continually increased. The most familiar SNIa datasets include, for instance, Gold04 Gold06 01 , 
SNLSij, ESSENCE [9l, SDSS [13, Davis [l^, and, most recently. Union Constitution 

In 414 SNIa from some heterogeneous compilations have been analyzed with the same analysis 
procedure. After selection cuts, this compilation was reduced to 307 SNIa, which have been named the 
Union dataset [ll|. There are 250 high redshift SNIa (z > 0.2) and 57 low redshift SNIa (z < 0.2) in the 
Union dataset ll|. Very recently, the CfA3 sample |13|] has been added to the 307 SNIa Union dataset 
to form the Constitution dataset [l^- Originally, the CfA3 sample consisted of 185 SNIa, which are all 
at the fairly low redshift, i.e., z < 0.08. After applying the same Union cuts, 90 SNIa from CfA3 survive. 
This increases the low redshift sample (z < 0.2) to 147 SNIa in the resulting Constitution dataset ^12l |. 
The 397 SNIa Constitution dataset is the largest published, spectroscopically confirmed sample to date. 

Soon after the release of the 397 SNIa Constitution dataset, some authors have used it to study dark 
energy. For example, the authors of ,14] smoothly reconstructed the deceleration parameter q{z) and the 
Om{z) diagnostic by using the Constitution dataset, and found that the cosmic acceleration might be 
slowing down. This is an unexpected result in some sense. The authors of ^15] compared the hologra phi c 
dark energy model , the Ricci dark energy model [iTj and the new agegraphic dark energy model jl8| , 
by using the Constitution dataset and other observational data; they found that the holographic dark 
energy model is more favored. This result is also somewhat different from the previous studies. Further, 
the authors of [l^ found that dark energy seemingly did not exist in the past and suddenly emerged at 
redshift z ^ 0.331, by fitting to the Constitution dataset alone. Needless to say, this is a striking claim. 
In addition, the authors of [20] found that the dark energy equation of state parameter (EoS) deviates 
from the cosmological constant at z > 0.5 significantly, by using the SNIa dataset (Union or Constitution) 
together with the gamma-ray bursts (GRBs) data; they noted that the result from Constitution dataset 
is somewhat different from others. In [2^, they also suggested that the deviation might arise from some 
biasing systematic errors in the SNIa and/or GRBs datasets. 

These unusual results bring our attention from dark energy to the SNIa dataset itself, especially the 
Constitution SNIa dataset. In fact, this situation is reminiscent of the one of Gold04 and Gold06 SNIa 
datasets, which also bring some interesting results. In [U [l^, the Gold04 dataset has been shown to 
be in 2a tension with the SNLS dataset and the WMAP observations. Although the Gold04 sample 
updated to Gold06 several years later, the tension still persisted. In [l^, the tension and systematics in 
the Gold06 SNIa dataset has been investigated in great detail. In the present work, we will follow the 
method used in ^] to study the tension in the recent SNIa datasets. 



II. OBSERVATIONAL DATA 



Before plunging into the issue of the tension in the recent SNIa datasets, here we briefly present the 
cosmological model and the observational data. As is well known, the SNIa data points are given in terms 
of the distance modulus i^obsizi). The theoretical distance modulus is defined as 

^J,th{Zl) = 5\ogiQDL{zi) + fio, (1) 

where no = 42.38 — 51ogio h and h is the Hubble constant Hq in units of 100 km/s/Mpc, whereas 

in which E = H/Hq and H is the Hubble parameter; p denotes the model parameters. In the present 
work, we consider the familiar Chevallier-Polarski-Linder (CPL) model [IJ, in which the EoS of dark 
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energy is parameterized as 



Wde ^ Wo + Wail ~ a) ^ Wo + Wa 



1 + z 



(3) 



where wo and Wa are constants. As is well known, the corresponding E(z) is given by [25l-l2 



E(z) 



+ + (1 - a^o) (1 + z)3(i+-o+».) 



exp 



ZWaZ 

l + z 



1/2 



(4) 



where Vl,no is the present fractional energy density of pressureless matter. The from the SNIa data is 
given by 



(p) = E 



[tJ'obsjzi) - Hth{zi)Y 

a\z,) 



(5) 



where tr is the corresponding Icr error. The parameter /iq is a nuisance parameter but it is independent of 
the data points. One can perform an uniform marginalization over /iQ. However, there is an alternative 
way. Following [2^ [28„ 2^ , the minimization with respect to no can be made by expanding the xfi of 
Eq. (O with respect to /zq as 



(6) 



where 



[^J'obsizi) - ^lth{zi\^J.o = o, p)] 



Eq. (|n]) has a minimum for j^io — B/C at 

X^(p)=i(p)-^. (7) 

Since xfi,min = xfi,min obviously, we can instead minimize which is independent of ^q. Note that the 
above summations are over the whole SNIa dataset. 

There are some other observational data relevant to this work, such as the observations of the cosmic 
microwave background (CMB) anisotropy 0, 01 and the large-scale structure (LSS) However, using 
the full data of the CMB and the LSS to perform a global fitting consumes a large amount of time. As an 
alternative, one can instead use the shift parameter R from the CMB, and the distance parameter A of 
the measurement of the baryon acoustic oscillation (BAO) peak in the distribution of SDSS luminous red 
galaxies. In the literature, the shift parameter R and the distance parameter A have been used extensively. 
It is argued that they are model- independent [soj , whereas R and A contain the main information of the 
observations of the CMB and the BAO, respectively. The shift parameter R is defined by [30,. .Slj 

where the redshift of recombination is z» = 1090, which has been updated in WMAP5 The shift 
parameter R relates the angular diameter distance to the last scattering surface, the comoving size of the 
sound horizon at z^ and the angular scale of the first acoustic peak in the CMB power spectrum of the 
temperature fiuctuations iSCj, IMI. The value of R has been updated to 1.710 ± 0.019 from WMAP5 
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The from the shift parameter R is Xb, — {R^ Robs)^ / <^r- On the other hand, the distance parameter 
A is given by 

where Zb = 0.35. In [11], the value of A has been determined to be 0.469 (ns/O.gS)-"-^'^ ± 0.017. Here the 
scalar spectral index Us is taken to be 0.960 which has been updated from WMAP5 3- The from the 
distance parameter A is = — AobsY /'^\- 

The best-fit model parameters are determined by minimizing the corresponding x^ • As in [25l . 133^ , the 
68% confidence level (C.L.) is determined by Ax^ = x^ — Xmm — l-^j 2.3 and 3.53 for = 1, 2 and 3, 
respectively, where Up is the number of free model parameters. Similarly, the 95% C.L. is determined by 
Ax^ = X^ - Xmvn < 4.0, 6.17 and 8.02 for 1,2 and 3, respectively. 



1 dS 
Zb Jo E{z) 



(9) 




FIG. 1: The 68% and 95% C.L. contours in the wq — Wa plane and the wq — QmO plane for the observations 
of SNIa only (blue dashed lines) and SNIa+yl+-R (black solid lines). We also show the best-fit values for the 
observations of SNIa only (blue box), SNIa-f-A (red star), SNIa+-R (green triangle) and SNIa+y4+_R (black point). 
These results are for the case of Constitution dataset. 



III. TENSION IN THE RECENT SNIA DATASETS 



At first, we consider the 397 SNIa Constitution dataset. We fit the CPL model to the observations of 
SNIa only, SNIa-|-A, SNIa-|-i?, and SNIa-l-^-fi?, respectively. The best-fit values are presented in Table HI 
In Fig. [1] we also present the 68% and 95% C.L. contours in the wq — Wa plane and the wq — plane. 
From Fig. [U we find that the Constitution SNIa dataset is in tension (significantly beyond 2(t) with the 
observations of the CMB and the BAO. 

Next, we turn to the 307 SNIa Union dataset. Similarly, we present the results in Table [TTl and Fig. [51 
From Fig. [2l we can see that the Union dataset is also in tension (beyond 2a) with the observations 
of the CMB and the BAO, although the corresponding tension is weaker than the one in the case of 
Constitution dataset. 
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Observation 


Xmin 




wo 




SNIa 


461.254 


0.453 


-0.207 


-11.316 


SNIa+yl 


465.438 


0.286 


-0.894 


-0.628 


SNIa+i? 


465.739 


0.260 


-0.930 


0.014 


SNIa+yl+i? 


466.1 


0.274 


-1.021 


0.413 



TABLE I: The Xmin B,nd the best-fit values of QmOi wo and Wa for the various observations. These results are for 
the case of Constitution dataset. 



Observation 






wo 


Wa 


SNIa 


310.091 


0.451 


-1.013 


-5.898 


SNIa+A 


310.915 


0.271 


-1.228 


1.502 


sm&+R 


310.911 


0.302 


-1.271 


1.269 


SNIa+A+7? 


311.154 


0.278 


-1.140 


0.859 



TABLE II: The same as in Table [Jl except for the case of Union dataset. 



Comparing the cases of the Constitution and Union datasets, we find that they have two common 
features: (i) the best-fit value of flmo from SNIa data only is fairly larger than 0.3; (ii) the best-fit value 
of Wa from SNIa data only is strongly negative, namely Wa is much smaller than —1. On the other hand, 
they have also two differences: (i) the Xmin 460 is fairly larger than the corresponding degrees of 
freedom dof ~ 397 for the Constitution dataset, whereas the Xmin ~ 310 is approximately equal to the 
corresponding dof ~ 307 for the Union dataset; (ii) for the case of Union dataset, the best-fit value of 
Wo remains about —1 for all observations, whereas it is not for the case of Constitution dataset. 

Now, one might ask whether or not all SNIa datasets are in tension with the observations of CMB and 
BAO. A third recent SNIa dataset is the Davis sample Jl3], which consists of 192 SNIa. We present the 
corresponding results in Table IIIII and Fig. |31 From Fig. |31 we see that the Davis SNIa dataset is fully 
consistent with the observations of CMB and BAO. There is no tension in the case of Davis dataset. For 
the case of 192 SNIa Davis dataset, the corresponding Xmin/'^'^f ^ 1' ^^'^ '"'o remains about —1 for 
all observations, similar to the case of Union dataset. However, its best-fit value of flmo from SNIa data 
only is much closer to 0.3 than the cases of Union and Constitution datasets. 



Observation 






wo 


Wa 


SNIa 


195.343 


0.349 


-1.105 


-1.229 


SNIa-fA 


195.453 


0.273 


-1.113 


0.526 


SNIa-fJ? 


195.463 


0.264 


-1.115 


0.664 


SNla+A+R 


195.485 


0.270 


-1.155 


0.818 



TABLE III: The same as in Table ID except for the case of Davis dataset. 
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FIG. 2: The same as in Fig. [T] except for the case of Union dataset. 




FIG. 3: The same as in Fig. [T] except for the case of Davis dataset. 



So, the situation is clear now: not all SNIa datasets are in tension with the observations of CMB and 
BAO. Besides the Davis SNIa dataset shown above, it is well known that the SNLS SNIa dataset is also 
fully consistent with the observations of CMB and BAO [2l|, Therefore, the recent SNIa datasets 
Constitution and Union are in tension not only with the observations of CMB and BAO, but also with 
the other SNIa datasets such as Davis and SNLS. In fact, they are in the similar situation of the Gold04 
and Gold06 SNIa datasets. 
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UnionOut subset (21 SNIa) 

1992bs, 1995ac, 1999bm, 1997o, 2001hu, 1998ba, 04Pat, 05Red, 2002hr, 03D4au, 
04D3cp, OSDlfc, 03D4dy, OSDlco, bOlO, d033, g050, g055, k430, ml38, m226 
ConstitutionOut subset (34 SNIa) 
1992bs, 1992bp, 1995ac, 1999bm, 1996t, 1997o, 1995aq, 2001hu, 1998ba, 04Pat, 05Red, 
2002hr, 03D4au, 04D3gt, 04D3cp, 03D4at, 03Dlfc, 04D3co, 03D4dy, 04D3oe, 04Dlak, 
03D1CO, bOlO, d033, f076, g050, k430, ml38, m226, snOlcp, sn02hd, sn03ic, sn07ca, sn07R 



TABLE IV: The names of SNIa in the UnionOut and ConstitutionOut subsets. 

IV. OUTLIERS WITH RESPONSIBILITY FOR THE TENSION 

In this section, we elucidate sources of the tension in the Constitution and Union SNIa datasets. Since 
the Constitution dataset comes from the Union dataset by adding 90 low redshift SNIa from the CfA3 
sample, and the tension exists in both Constitution and Union datasets, we speculate that the main 
sources might be in the Union dataset, and the tension in the Constitution dataset is mainly inherited 
from the Union dataset. To verify this speculation, we divide the Constitution sample into the high 
and low redshift groups with the same dividing line z — 0.2 in [ill, [l^. and repeat the procedure in the 
previous section. We find that the results for the high and low redshift groups have only little difference, 
and this supports our speculation mentioned above. So, we firstly pay our attention to the Union dataset. 

In the case of Gold06 SNIa dataset, there are only five subsets (see Table I of [11]). Using the method 
of subset truncation (see Sec. IVl below) . the authors of 23;] found that the HZSST subset is the main 
source of the tension in the Gold06 dataset. Then, they isolated six SNIa in the HZSST subset which are 
mostly responsible for the tension. On the side of Union dataset, there are 13 subsets in it (see Table 3 
and Table 11 of [HI). Even throwing away the first five subsets at low redshift, there are still eight 
subsets to be considered. By careful observation, we find that three subsets (Barris, Perlmutter, and 
Riessl998+HZT) are notable, because their RMS listed in Table 3 of [IH are significantly higher than 
other subsets. Then we try to subtract one of these three subsets in turn, and subtract these three subsets 
together from the full Union dataset, and see whether or not the tension can be removed. However, we 
find that it does not work in fact. This leads us to speculate that the main sources of the tension in 
the Union dataset might be homogeneously distributed in the whole Union dataset, namely, they do not 
concentrate in a single subset listed in Table 3 of [ll|, unlike the case of Gold06 dataset. 

Here, we follow the simple method used in [2^ to find the outliers responsible for the tension. In [23j . 
the distance moduli of the six SNIa which are mostly responsible for the tension in Gold06 dataset differ 
by more than l.Scr from the ACDM (with fimo = 0.28) prediction. Similarly, we firstly fit the flat ACDM 
model to the whole 307 SNIa in the Union dataset, and find that the best-fit parameter is rimo = 0.287 
(the corresponding /io = BjC = 43.16). Then, we calculate the relative deviation to the best-fit ACDM 
prediction, ]/io6s ~ [i-KCDuM ^^ohs^ for all the 307 points. There are 16 SNIa which differ from the best- 
fit ACDM prediction beyond 2f7. However, we find that the tension cannot be completely removed by 
subtracting these 16 SNIa from the 307 SNIa Union dataset. We should adopt a severer cut. There are 
21 SNIa differing from the best-fit ACDM prediction beyond 1.9cr, and we call them the "UnionOut" 
subset. The names of these 21 SNIa are listed in Table ITVl In Fig. |4j we also plot their distance modulus 
deviations relative to the best- fit ACDM prediction. These 21 SNIa are indeed homogeneously distributed 
in the whole Union dataset. On the other hand, their minimum and maximum redshift are Zmm = 0.0488 
(1995ac) and z^ax — 1.1900 (05Red), respectively. It is worth noting that if one chooses to use instead 
the XCDM model in which w^e = const, to find the outliers, the results are almost the same of the case 
of ACDM model used here. 

By subtracting the 21 SNIa UnionOut subset from the whole 307 SNIa Union dataset, we obtain the 
so-called "UnionT" sample ("T" stands for "truncated"). Obviously, the UnionT sample consists of 286 
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UnionOut subset 



0.2 0.4 0.6 O.i 



1 1.2 




FIG. 4: A/i — fiobs(zi) — i-iACDM{zi) for the SNIa in the UnionOut and ConstitutionOut subsets, where the best-fit 
ACDM are for the whole Union and Constitution datasets, respectively. See the text for details. 



Observation 


Xmin 




Wo 




SNIa 


204.060 


0.360 


-1.250 


0.242 


SNIa+^ 


204.067 


0.277 


-1.151 


1.206 


SNIa+i? 


204.067 


0.298 


-1.181 


1.064 


Sma.+A+R 


204.194 


0.281 


-1.086 


0.735 



TABLE V: The same as in Table U except for the case of the 286 SNIa UnionT sample. 



SNIa. Then, we repeat the same fitting as in Sec. IIIII for this 286 SNIa UnionT sample, and present the 
corresponding results in Table |V] and Fig. [SJ From Fig. [51 we clearly see that the UnionT SNIa sample is 
fully consistent with the other observations, and the tension has been completely removed. On the other 
hand, by dropping the 21 outliers, the Xmin significantly reduces from Xmin ~ 310 for the 307 SNIa Union 
dataset (see Table HI)) to Xmin 204 for the 286 SNIa UnionT sample (see Table IV)) : the corresponding 
Xmin/'^'^f been improved. 

Notice that since dropping these 21 SNIa from the Union dataset is enough, we need not adopt other 
severer cuts, such as 1.8cr or 1.7a. To preserve the number of usable SNIa as much as possible, the cut 
1.9a adopted here is the best choice. 

Now, let us turn to the Constitution dataset. Similarly, we firstly fit the flat ACDM model to the 
whole 397 SNIa in the Constitution dataset, and find that the best-fit parameter is fimo = 0.290 (the 
corresponding fig — B/C = 43.3158). Then, we calculate the relative deviation to the best-fit ACDM 
prediction, \fiobs — ^J'ACDM\/o^obs, for all the 397 points. We adopt the same cut 1.9a used in the Union 
dataset. There are 34 SNIa differing from the best-fit ACDM prediction beyond 1.9a, and we call 
them the "ConstitutionOut" subset. The names of these 34 SNIa are listed in Table ITVl In Fig. SI we 
also plot their distance modulus deviations relative to the best-fit ACDM prediction. These 34 SNIa 
are indeed homogeneously distributed in the whole Constitution dataset. On the other hand, their 
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FIG. 5: The same as in Fig. \T\ except for the case of the 286 SNIa UnionT sample. 



minimum and maximum redshift are Zmin = 0.015 (sn07ca) and Zmax — 1.190 (05Red), respectively. 
Obviously, the UnionOut subset and ConstitutionOut subset do heavily overlap. In fact, most SNIa of 
the ConstitutionOut subset come from the Union sample. Only the last five SNIa in the ConstitutionOut 
subset come from the CfA3 sample. This confirms our speculation in the beginning of the present section 
that the tension in the Constitution dataset is mainly inherited from the Union dataset. In other words, 
the tension is not caused by the low redshift CfA3 sample. 

By subtracting the 34 SNIa ConstitutionOut subset from the whole 397 SNIa Constitution dataset, we 
obtain the so-called "ConstitutionT" sample. Obviously, the ConstitutionT sample consists of 363 SNIa. 
Then, we repeat the same fitting as in Sec. IIIII for this 363 SNIa ConstitutionT sample, and present 
the corresponding results in Table I VII and Fig. [HI From Fig. [SI we see that the ConstitutionT SNIa 
sample is consistent with the other observations at the 2cr level, and the tension has been significantly 
alleviated, especially comparing Fig. [6] with Fig. [T] On the other hand, by dropping the 34 outliers, 
the Xmin significantly reduces from Xmin ~ 461 for the 397 SNIa Constitution dataset (see Table |I| to 
Xmin ~ 269 for the 363 SNIa ConstitutionT sample (see Table |VI|; the corresponding Xmin/^'^f been 
significantly improved. 

Notice that if one instead adopt a severer cut 1.8a for the case of Constitution dataset, the results are 
close to the one of 1.9tT in fact. So, to preserve the number of usable SNIa as much as possible, the cut 
1.9(T adopted here is appropriate. 



Observation 


Xmin 




wo 


Wa 


SNIa 


268.9 


0.407 


-0.988 


-2.270 


SNIa+A 


269.13 


0.286 


-0.983 


0.383 


sm&+R 


269.134 


0.282 


-0.984 


0.438 


SNla+A+R 


269.138 


0.284 


-0.996 


0.486 



TABLE VI: The same as in Table jT] except for the case of the 363 SNIa ConstitutionT sample. 
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FIG. 6: The same as in Fig. [T] except for the case of the 363 SNIa ConstitutionT sample. 



V. VERIFYING THE OUTLIERS WITH THE METHOD OF RANDOM TRUNCATION 

In the previous section, we have identified the outhers which are mostly responsible for the tension in 
the Union and Constitution datasets. Here, we will make them more robust by employing the method of 
random truncation used in 23]. 

First, we consider the Union SNIa dataset. Following [2^, we compare the best-fit values of ilmo, wq 
and Wa for the 286 SNIa UnionT data with a large number A'^ of corresponding random truncations of the 
Union SNIa data (here we adopt the same number = 500 as in [1^). The random truncations involve 
random subtractions of the same number of SNIa and in the same redshift range as the UnionOut subset 
from the full Union dataset (notice that UnionT ~ Union — UnionOut). We can easily obtain the mean 
best-fit value and the la range aqr of the quantity q = Qmo, wq and Wa for these random truncations. 
According to |23], if the best- fit quantity q of the UnionT sample is within the la range of the mean 
best-fit value of the random truncations, the UnionOut subset is a typical truncation representative 
of the Union dataset and statistically consistent with it. If on the other hand the best-fit quantity q of 
the UnionT sample differs from the mean best-fit value q^ of the random truncations beyond 2a, we can 
conclude that the UnionOut subset is not a typical truncation and is systematically different from the full 











q-q'' 




= 0.360 




= 0.448 ± 0.025 


-3.6cr 


Wo = 


-1.250 


w5 = 


-1.028 ±0.153 


-1.4cr 


Wa = 


= 0.242 




-5.715 ± 1.902 


4-3.1(7 



TABLE VII: The first column is the best-fit values of the quantity q for the 286 UnionT SNIa sample, which can 
be read from Table |Vl The second column is the mean best-fit value and corresponding la range of the quantity 
q for the A'^ = 500 random truncations. The third column is the relative deviation. See the text for details. 



11 



& -5 



-10 



-15 
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Wo 
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0.5 



0.4 



0.3 



0.2 



-2 -1.5 -1 -0.5 0.5 

Wo 



FIG. 7: Comparing the best-fit parameters of the 286 UnionT SNIa sample (red stars) with the ones of the 500 
random truncations (orange points) in the mo — Wa plane and the wq — S^mo plane. The mean values of these 500 
points are also indicated by the black solid points. See the text for details. 



Union dataset [1^ . We present the results in Table IVIII and Fig. [T] Obviously, the UnionOut subset are 
systematically different from the full Union dataset, and its 21 SNIa are indeed the outliers. Therefore, 
it is reasonable to cut these 21 SNIa outliers in the UnionOut subset from the full Union dataset. 

Next, we turn to the case of the Constitution dataset. Similarly, we follow the procedure described above 
(notice that ConstitutionT = Constitution — ConstitutionOut), and present the results in Table Ivnll and 
Fig. [51 Significantly, the ConstitutionOut subset are systematically different from the full Constitution 
dataset, and its 34 SNIa are indeed the outliers. Therefore, it is reasonable to cut these 34 SNIa outliers 
in the ConstitutionOut subset from the full Constitution dataset. 



VI. CONCLUDING REMARKS 



In the present work, we investigated the tension in the recent SNIa datasets Constitution and Union. 
We have shown that they are in tension not only with the observations of CMB and BAO, but also 
with other SNIa datasets such as Davis and SNLS. Then, we found the main sources which are mostly 
responsible for the tension. Further, we made this more robust by employing the method of random 
truncation. 

Based on the results of this work, we suggest to perform the severer selection cut in the Union and 
Constitution SNIa datasets to reject the outliers further. While one uses the full 307 SNIa Union dataset 
and the full 397 SNIa Constitution dataset to constrain the cosmological models, we strongly recommend 
use also of the 286 SNIa UnionT sample and the 363 SNIa ConstitutionT sample given in this work for 
comparison. In fact, it is anticipated that the results might be fairly different, and the unusual features 
which arise from the full 307 SNIa Union dataset and the full 397 SNIa Constitution dataset might 
disappear in the cases of the 286 SNIa UnionT sample and the 363 SNIa ConstitutionT sample. 

As communicated by [1^, the era of precision cosmology has arrived (e.g. [34] and [Tl) and there is a 
deemed need for an increased understanding of fundamental supernova (SN) physics. In addition, future 
SN datasets, e.g., those from DES and LSST, will rely significantly on photometric redshift determinations 
and SN type classifications (particularly in the case of the LSST). The details of SN colors will have to be 



12 



= 0.407 
Wo = -0.988 
Wa = -2.270 



= 0.453 ± 0.011 
= -0.198 ±0.171 
w': = -11.470 ± 1.804 



-4.0(7 
-4.6(7 

+ 5.1(7 



TABLE VIII: The first column is the best-fit values of the quantity q for the 363 ConstitutionT SNIa sample, 
which can be read from Table IVTl The second column is the mean best-fit value and corresponding la range of 
the quantity q for the A'' — 500 random truncations. The third column is the relative deviation. See the text for 
details. 
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FIG. 8; Comparing the best-fit parameters of the 363 ConstitutionT SNIa sample (red stars) with the ones of the 
500 random truncations (orange points) in the wo — Wa plane and the wo — ^mo plane. The mean values of these 
500 points are also indicated by the black solid points. See the text for details. 

unraveled if purely photometric SN are to be successfully used for precision cosmology [13] . Enlightened 
by the results of this work, one might further consider the possible impact of the ConstitutionOut and 
UnionOut subsets on the understanding of fundamental SN physics and SN colors [35i|. We leave this 
issue as an open question. 
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Note added: ^fter the submission of the present paper, the 1st year SDSS-II supernova sample has 
been released |34|. It consists of 103 SNIa with redshifts 0.04 < z < 0.42. Here, we would like to have 
a primary analysis. Note that the 103 SDSS SNIa dataset was given in both the frameworks of using 
SALT2 and MLCS2k2 light curve fitters. Similarly, we fit the CPL model to the observations of SNIa 
only, SNIa+j4, SNIa+i?, and SNIa+A+i?, respectively. The best-fit values are presented in Table IIXI 
We can see that the best-fit parameters flmo, wq and Wa for the case of SNIa only (both SALT2 and 
MLCS2k2) are fairly unusual. On the other hand, we find that in the case of SNIa only, the constraints 
on the model parameters flmo, wq and Wa are very loose. These results might be mainly due to the 
relatively narrow redshift range (0.04 < z < 0.42) and the relatively small number of the SDSS SNIa 
sample. Considering these issues, unlike the cases of the Constitution, Union and Davis datasets, the 
tension analysis might be not so robust for the case of the SDSS SNIa dataset. However, since the SDSS 
SNIa dataset fills in the redshift "desert" between low- and high-redshift SN surveys, it deserves further 
attentions in the SN studies. 



Observation 


Xmin 




■wo 




SNIa (SALT2) 


122.305 


0.999 


90.379 


-1776.16 


SNIa (SALT2) + A 


128.902 


0.297 


-1.081 


2.147 


SNIa (SALT2) + R 


128.909 


0.406 


-1.234 


1.417 


SNIa (SALT2) + A + R 


128.973 


0.296 


-0.881 


0.252 


SNIa (MLCS2k2) 


173.159 


0.970 


29.806 


-438.023 


SNIa (MLCS2k2) + A 


181.143 


0.315 


1.497 


-19.822 


SNIa (MLCS2k2) + R 


181.646 


0.137 


1.053 


-14.190 


SNIa (MLCS2k2) + A + R 


194.408 


0.311 


-0.532 


-0.877 



TABLE IX: The same as in Table HI except for the case of SDSS SNIa dataset. 
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